Collision rates in near-resonant optical lattices 
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We present a simple method to calculate the binary collision rate between 
atoms in near-resonant optical lattices. The method is based on the Monte 
Carlo wave function simulations and the collision rate is obtained by monitoring 
the quantum flux beyond the average distance between the atoms. To illustrate 
the usefulness of the method, we calculate the collision rates for a wide range 
of occupation densities and various modulation depths of the lattice. The 
method presented here combined with the semiclassical calculations accounting 
for intra- well collisions can simplify the study of the effects of binary collisions 
on the dynamics of atomic clouds trapped in near-resonant optical lattices. 
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1. Introduction 

A periodic polarization or intensity gradient of a laser 
field can create a periodic optical potential structure for 
atoms. When atoms move in this structure, they may 
undergo Sisyphus 'rftJ^Snfpd finally localize into the 
optical lattice sitesEMaafflll. After localization, atoms 
are still able to move around in the lattice. The optical- 
pumping-induced motion typically dominates in near- 
resonant optical lattices, whereas in shallow far-off reso- 
nant lattices the quantum mechanical tunneling of atoms 
between the lattice sites may dominate the atomic mo- 
tion. Because of their ability to move between the lattice 
sites, two atoms may end up in the same site and collide. 
The purpose of this paper is to present a simple method 
to calculate the binary collision rate in a near-resonant 
optical lattice. 

Cold collisions have been widely studied in magneto- 
optical trapaj'EI. This is not the case for optical lattices 
where the complications arise due to the position depen- 
dent coupling between the multi-level atoms and the laser 
field. So far, to the best of our knowledge, there has 
been only a few., experimental cold collision studies in 
optical latticesE3Eil. Theoretical studies of the interac- 
tions between atoms in.pjaticaj. lattices include the mean- 
field type approachesOE3EJIljll3, and Monte Carlo wave 
function (MCWF) simulations of binarji-collisions in red 
detunecOEa and blue detuned latticesEj. However, the 
mean-field approaches neglect the dynamical nature of 
the collisions, and the MCWF simulations for two atoms 
colliding in optical lattice require extremely large com- 
puter resourcealj. It would therefore be useful to find 
a way to do a dynamical study of collisions in near res- 
onant optical lattices without the requirement of heavy 
computational resources. 

We present here a simple method to calculate the 
rate of binary collisions in near-resonant optical lattices. 
The .mpthad .is ..based on single- atom MCWF simula- 
tionsEJOESSoLJ and the key idea is to monitor a quan- 
tum flux beyond the average distance z a between the 



atoms in a densely populated lattice. The accumulation 
of the atomic population beyond z a gives information 
about the atomic quantum transport in a lattice and can 
be monitored by MCWF simulations for one atom. Thus, 
our method avoids the requirement of large computa- 
tional resources of the two-atom MCWF collision simu- 
lationMM 

To illustrate the usefulness of the method we calcu- 
late the collision rates for various modulation depths and 
occupation densities of one dimensional optical lattices. 
Moreover, we show that it is possible to obtain accurate 
collision rates for all the densities from only a few sim- 
ulations for a given lattice. The simulation results show 
a quadratic behaviour of the binary collision rate as a 
function of density. 

It has been shown before that, in the parameter regime 
we use here, the atomic motion between the lattice sites 
is to a good anproximation dominated by the laser-atom 
interactions Ota. The effects of the interactions between 
atoms come into play as soon as the two atoms try to oc- 
cupy the same lattice site and collide. Thus the combina- 
tion of the method presented here to calculate the binary 
collision rate, and the semiclassica.J-|tnethods accounting 
the effects of intra- well collisions £3oc3, can simplify the 
study of the effects of collisions on the dynamics of the 
atomic cloud trapped in a near-resonant optical lattice. 

We present the lattice structure and our calculation 
method of the collision rate in the next section, the re- 
sults for various lattice depths and occupation densities 
in Section ||, and conclude with the discussion in Section 

§ 

2. Optical lattice and collision rate 

A. Sisyphus cooling in optical lattice 

We consider here atoms having ground state angular mo- 
mentum J g = 1/2 and excited state angular momentum 
J e = 3/2 corresponding to alkali metal elements when the 
hyperfine structure is neglected, and use the atomic mass 
M of 133 Cs. The resonance frequency between the states 
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Table 1. Used parameters. Rabi frequency fl, detuning 
5, lattice modulation depth Uq, and saturation parameter 

«o- 

fi(r) 8(F) U {E r ) s 



1.0 
1.5 
1.9 



-3.0 
-3.0 
-3.0 



259 
580 
936 



0.054 
0.122 
0.195 



is ujq so that Hloq = E e — E g , where E e and E g are ener- 
gies of the ground and the excited states in zero field. A 
single atom has two ground state sublevels \g±±/2 > and 
four excited state sublevels |e±3/2 > and \e±i/ 2 >, where 
the half-integer subscripts indicate the quantum number 
m of the angular momentum along the z direction. 

The laser field consists of two counter-propagating 
beams with orthogonal linear polarizations and with fre- 
quency u>l. The total field has a polarization gradient in 
one dimension and reads (after suitable choices of phases 
of the beams and origin of the coordinate system) 



E(z,t) = £ (e a 



Ak r z 



-ik, 



z )e 



-iuj L t 



+ C.C. 



(1) 



where £$ is the amplitude and k r the wavenumber. 

The intensity of the laser field and the strength of the 
atom-field coupling is described by the Rabi frequency 
51 = 2d£o/h, where d is the atomic dipole moment of 
the strongest transition between the ground and excited 
states. The detuning of the laser field from the atomic 
resonance is given by 5 = lol — uiq. As a unit for fl and 
S we use the atomic linewidth V, and express energy in 
the recoil unit E r = (h 2 k 2 )/2M. 

We keep the detuning fixed, S = — 3T, and vary the 
Rabi frequency Q, which gives various values for the op- 
tical potential modulation depth 



U 



where sq is the saturation parameter given by 

n 2 /2 



so 



S 2 + T 2 /A' 



(2) 



(3) 



See Table |l] for used parameters. 

The system Hamiltonian after the rotating wave ap- 
proximation reads 



H s = ^—- hSP e + V. 
2M 



(4) 



is the 

n) (6m|j 

atom and the field is 



Here 



Pe 

-3/2 



V = 



-i — — sim kz ) 
y/2 

— ■ = cos(fez) 
V2 

-h.c. 



projector operator P e = 
and the interaction between a 



|e 3 /2> (51/2I + -7=|ei/ 2 ) (ff-i/2| 



-3/2) 



V3> 

(5-1/2I + ^|e_i/ 2 ) (9 



-1/2 



The position and time dependent wave function of the 
system is 



\i>{z,t)) =^2<i/) jtm (z,t)\jm) 



(6) 



The periodic polarization gradient of the laser field and 
the subsequent periodicity of the couplings between the 
atomic states are reflected in the periodic light shifts (ac- 
Stark shifts) of the atomic sublevels creating the optical 
lattice structure. When the atomic motion occurs in a 
suitable velocity range, optical pumping of the atom be- 
tween the groundpStatei-sublevels reduces the kinetic en- 
ergy of the atom tyflfi'QEBQ, and the atom is cooled. 

When the steady state is reached after a certain period 
of cooling, atoms are to a large extent localized into the 
optical potential wells. In this study we deal with near- 
resonant optical lattices where the laser field is detuned a 
few atomic linewidths to the red of the atomic transition. 
Thus, after localization atoms may still move around in 
the lattice due to the finite extent of the atomic wave 
packet and scattering of photons which transfers the pop- 
ulation between the various ground state sublevels and 
corresponding optical potentials via the optical pumping 
mechanism. 

Because the atoms are able to move between the lattice 
sites, they may try to occupy the same site and collide. 
The purpose of this paper is to find a simple way to calcu- 
late the binary collision rate without the need to perform 
complicated two-atom simulations. We use the occupa- 
tion densities po of the lattice between 14.2% and 25.0% 
corresponding to cases with every seventh and fourth site 
occupied respectively. Other relevant parameters are pre- 
sented in Table [l]. From the simulation results, it is pos- 
sible to calculate the collision rates for a wide range of 
occupation densities of the lattice, see Section [|. 

B. Calculation of collision rate 
1. Monitoring the quantum flux 

The key idea in our method is to monitor the quantum 
flux beyond the average distance z a between the atoms in 
the one dimensional lattice. Since the distance between 
the nearest neighbour lattice sites is A/4, then, e.g., z a — 
A corresponds to the occupation density po — 25% of the 
lattice. Here A is the wavelength of the lattice lasers. 
We call the area \z\ > z a the accumulation region (see 
Fig. |]) and denote the time dependent fraction of the 
wave packet in this region by 



m 2 (t) = 



ip*(z, t)il)(z, t)dz. 



(7) 



(5) 



The ensemble averaged Monte Carlo (MC) result of 
|?/> a | 2 (£) gives information about the number of collisions 
and the collision rate, as shown below. 

To get the value for a collision rate we need to know 
how many atoms have travelled into the accumulation 
region and in which time. Thus, we need information 
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about the quantum flux into the accumulation region, in 
other words the cumulative population in the accumula- 
tion region as a function of time. This means that we 
should be able to make the quantum flux unidirectional 
into \z\ > z a . The atoms should therefore be able to 
arrive into the accumulation region due to their random 
motion in the lattice but should not be allowed to leave. 
In this case the accumulation rate of the atomic popula- 
tion into the accumulation region could give us the binary 
collision rate in the lattice. 

To make the quantum flux unidirectional into the ac- 
cumulation region we add the diagonal term 

0:|z| <z a -A/8, 
—a :\z\ > z a + A/8, 
• + 1/8)] :-z a - A/8 < z ^ 

-z a + A/8, W 
■ + 1/8)] :z a - A/8 s$ z ^ 
z a + A/8 



H„ 



—a cos 2 [2tt( 
-a sin 2 [2n( 



. -\- ; , 
A 



A 



to the system Hamiltonian, Eq.([|). We show the total 
optical potentials for the two ground state sublevels in 
Fig. |l| The steady state momentum distributions and 
the average kinetic energy of the atoms in the region 
\z\ < z a agree with the already known results for optical 
lattices, see for example Ref. 27. 

When an atom arrives into a lattice site located around 
±z Q , corresponding to a collision site, it is pushed into 
the accumulation region and is not allowed to climb back 
to the region \z\ < z a . This does not change the physics 
for region \z\ < z a before the atom reaches collision site. 
After arrival the atom is captured into the accumulation 
region and as a consequence it is easy to calculate the 
cumulative population as a function of time. The trick 
here is that the drop in energy can not be too large but 
has to be large enough: Sisyphus cooling should still be 
effective in accumulation region, preventing the atoms 
from simply bouncing back from the wall of the grid and 
traveling back to \z\ < z a \ but the energy drop has to be 
large enough so that the probability to climb the energy 
barrier from the accumulation region should be small. 
We have chosen the values a = 375, 840, 1356EV for Uo — 
259, 580, 936.B,. respectively. 

The total Hamiltonian for the time evolution of the 
system now reads 



H — H s + H n 



Ho- 



rn 



where H s is the system Hamiltonian from Eg. (Q), H m 
the optical potential modification from Eq. (^), and Hd 
includes the non-Hcrmitian decay part of the MC method 
(see below). 

2. Monte Carlo wave function method 

We use the MCWF method (see Refs. 20-23. for de- 
tails, and, e.g., Refs. 24-26 for applications of the method 
to the cold collision problems) to calculate the collision 
rates. We have also applied the method recently to study 
heating in red-detuneorlri and optical shielding in blue- 
detuned latticestUl 




Fig. 1. The schematic view of the optical potentials 
for the two ground state Zeeman sublevels. The initial 
lattice site of the wave packet is indicated by the arrow 
and the modification of the potentials around the points 
±z a makes the quantum flux unidirectional into the ac- 
cumulation region (shown as shaded area). 



The core idea of the MCWF method is to generate 
a large number of single wave function time evolution 
histories which include stochastic quantum jumps. The 
jumps occur because of the non-Hermitian part Hd in the 
system Hamiltonian. The information about the decay 
of the system is included in Hd, which shrinks the norm 
of the wave functionj-and gives the jump probability for 
each time step takenca. The results for the system prop- 
erties are finally obtained as an ensemble average of the 
generated single wave function histories. 

The initial position of the wave packet in our simula- 
tions is given randomly into the lattice well around the 
point z = (see Fig. [fl), and with zero mean momen- 
tum. Strictly speaking we should use as initial position 
and momentum distributions the steady state distribu- 
tions for the used lattice parameters. We emphasize that 
the time scale to achieve the steady state is short com- 
pared to the population accumulation time to the accu- 
mulation region and we have checked that the change 
of the initial conditions does not change the simulation 
results. Thus, we avoid doing the double effort (first cal- 
culating the steady state properties of the system and 
then the collision rate) but still get the correct results. 
If one does simulations for higher occupation densities 
than used here (which is not actually necessary since the 
simulation results calculated here can be used to obtain 
the rate curve for wide range of densities, see Section ||), 
then more care should be given to the initial conditions 
of the system. 

The number of the generated wave function histories 
per simulation varies from 256 to 320 and the total sim- 
ulation times vary from 3200r _1 to 6400r~ 1 . 
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Fig. 2. An example (U = 936E r:Po = 20%) of the ac- 
cumulation region population \ipa\ 2 as a function of time. 
The binary collision rate R is obtained from the slope of 
the curve f3 and the average distance between the atoms 
z a as R — P/z a , see text. 




p (%) 



Fig. 3. The binary collision rate R for three different 
lattice depths Uq as a function of occupation density po 
of the lattice. The points show the simulation results and 
the solid lines the quadratic collision rate curves averaged 
from the simulation results for the specific lattice. 



3. Collision rate 

We show an example of the ensemble averaged atomic 
population in the accumulation region, | -0 a | 2 , as a func- 
tion of time in Fig. || where we display the simulation 
result and the linear fit. 

Figure || demonstrates a steady flow of atomic popu- 
lation into the accumulation region. Thus, by fitting a 
linear function to the MC simulation result we can use 
the slope of the curve, (3, to calculate the collision rate 
R in the lattice. 

In the example of Fig. | \ip a \ 2 {t = 3000F- 1 ) » 0.4 
with the number of histories of 320. This means that the 
number of collisions in time t — 3000r _1 in this specific 
case was 0.4 per z a or 0.4 x 320 = 128 for a lattice length 
of 320z a . 

In general, the collision rate from simulations per unit 
time and per unit volume (per unit length in our one 
dimensional lattice) is given by 

R = p/z a . (10) 

The total number of collisions N to t hi the experiment 
would thus be N to t — R x t x L where t is time and L 
the length of the one dimensional lattice. 

3. Results 

We have simulated the binary collision rate R for three 
different lattice depths Uq with fixed detuning S = 
— 3r (see Table [l]), and for lattice occupation densities 
po = 14.3%, 16.7%, 20%, 25%, from every seventh to ev- 
ery fourth lattice site occupied respectively. The results 
are displayed in Fig. ||. 

The simulation results showr&,quadratic behaviour of R 
with respect to atomic densityES. Thus we can obtain the 
whole collision rate curve for all densities of the specific 



lattice by calculating R from only a single density MC 
result. By dividing the result for R with the square of the 
density it is possible get the factor of quadratic behaviour 
for all densities. We note that this might be vulnerable to 
statistical error typical for MC simulations. To improve 
the statistical accuracy of the result for all densities we 
take the average of the quadratic factors obtained for 
the different densities. This makes it possible to obtain 
the whole collision rate curve for all the densities of the 
specific lattice with good statistical accuracy. The final 
results for all the density range obtained by averaging 
the various density MC results are drawn with solid lines 
in Fig. |. 

In the parameter regime we use, the motion between 
the lattice sites is mainly due to optical pumping, not 
because of the tunneling between the lattice sites. The 
higher is Uq, the larger is the collision rate. For fixed de- 
tuning the increasing field intensity means higher f7 and 
shorter optical pumping time. The atoms get more mo- 
bile since the increasing excited state population makes 
the probability of the internal state changing sponta- 
neous emission event higher. If the atomic motion be- 
tween lattice sites was because of tunneling, then the 
collision rate would decrease for tightly bound atoms for 
larger lattice modulation depths. 

Moreover, the atomic motion between the lattice sites 
is not simply a random walk with a step size of the lat- 
tice constant. When the atom changes its internal state 
and optical potential well, it is not necessary that the 
atom halts its motion already in the nearest neighbour 
site. The atom may travel the distance of several lattice 
constants with a single flight, in fact, for the shallow opti- 
cal lattices the atomic motion and diffusion in the lattice 
may even be anomalous and include Levy walksJIaLjiJ. 
For the parameters we use, the diffusion is norma£3c3. 
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Since the atomic motion between the lattice sites here 
is due to optical pumping it is useful to compare the 
ratios of calculated collision rates of the various lattice 
depths to the corresponding ratios of the optical pump- 
ing times. The ratio of the optical pumping times for 
two different lattice depths with equal detuning is sim- 
ply given by the square of the ratio of Rabi frequenciesEZI. 
Here the ratio of the optical pumping times is on the 
same order of magnitude as the calculated ratios of the 
collision rates, but there can be a difference of factor of 
two. They do not match exactly because there is more 
to atomic motion between the sites than optical pump- 
ing only. For example the hopping statistics between the 
sites may differ for various lattice depths. 

4. Discussion and conclusions 

We have presented a method to calculate the binary colli- 
sion rate R in a steady state for an atomic cloud trapped 
in a near-resonant optical lattice. The method is based 
on MCWF simulations and the key idea is to monitor the 
quantum flux beyond the average distance between the 
atoms. 

We have done the MCWF simulations for near- 
resonant lattice, 5 = —3T, and for lattice depths 259E r ^ 
U ^ 936E r . From the MCWF results it is possible to 
calculate the quadratic collision rate curve for all the den- 
sities of the lattice, within some limitations though (see 
below) . 

The advantage of the method is the avoiding of the 
large computational resource requirements of two-atom 
simulations. This is possible because in near-resonant 
lattices with large enough modulation depths the atomic 
motion between lattice sites is dictated by the field-atom 
interactions only EZHIil The interactions between the 
atoms come into play only when two atoms try to oc- 
cupy the same lattice site and collide (this is not the case 
for all parameter ranges, see Ref. 13). The method is also 
fast, straightforward, and fully quantum-mechanical. 

In addition of giving quantitative estimate for the num- 
ber of collisions in the experiment for a wide density 
range of optical lattices, our method brings out the pos- 
sibility to combine two simple and computationally light 
methods to study collisions in optical lattices. Namely, 
the method presented here, combined with the semiclas- 
sical c^Lcailations accounting for the intra-well collision 
effectsffle3, would simplify the study of the effects of the 
dynamic binary collisions for the atomic clouds trapped 
in optical lattices. 

We have actually calculated here the collision rate for 
the primary collisions in the lattice. That is: the atomic 
cloud achieves the steady state and it is the steady state 
properties which define the collision rate here. Naturally 
the dynamics of the cloud may change when a large frac- 
tion of the atoms collide in short time for high occupation 
densities of the lattice. Moreover, the occupation density 
may change if the collided atoms gain enough kinetic en- 
ergy to leave the lattice, or if the lattice is constructed 



for metastable atomstStLlI, which may ionize and escape 
the lattice when colliding. This is actually a benefit when 
using metastables, since it would be straightforward, at 
least in principle, to compare the number of collisions 
given by the method presented here and in the experi- 
ment. 

We have to neglect the rescattering of photons. It is 
practically impossible account the rescattering effects be- 
cause of finite amount of available computer resources. 
Despite of this, we think that our study is useful because 
of two reasons: a) we can give results for collision rates in 
the density region where the rescattering effects do not 
appear, b) because for very high densities (where rescat- 
tering effects should be accounted in principle) one could 
still use the method presented here to study the aspects 
of the binary collision effects on the lattice dynamics of 
the atomic cloud. This l b j i ts olf is very complex problem, 
and also worth studyingl3ll3't3. 

From a Monte Carlo method point of view there is an 
interesting feature present here. Namely, we first do a 
few simulations and notice that the binary collision rate 
has a quadratic behaviour with respect to the atomic 
density. As a second step, we obtain the final results by 
taking the average of the various sets of MC results. We 
emphasize that this method differs from a simple increase 
of the number of histories in the MC ensemble to increase 
the statistical accuracy. Thus the benefits here are (in 
addition to the ones mentioned in previous paragraphs) 
twofold a) the statistical accuracy of the results increases 
b) the result can be obtained in a wider range than in 
which the MC simulations are done. This is a new feature 
in MCWF simulations to our knowledge, at least when 
MCWF method is applied to cold collision problems. 
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